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Abstract 

We compute at strong coupling the large A^ correlation functions of supersymmet- 
ric Wilson loops in large representations of the gauge group with local operators 
of AA = 4 super Yang- Mills. The gauge theory computation of these correlators 
is performed using matrix model techniques. We show that the strong coupling 
correlator of the Wilson loop with the stress tensor computed using the matrix 
model exactly matches the semiclassical computation of the correlator of the 
't Hooft loop with the stress tensor, providing a non-trivial quantitative test of 
electric-magnetic duality of AA = 4 super Yang- Mills. We then perform these 
calculations using the dual bulk gravitational picture, where the Wilson loop is 
described by a "bubbling" geometry. By applying holographic methods to these 
backgrounds we calculate the Wilson loop correlation functions, finding perfect 
agreement with our gauge theory results. 
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1 Introduction and outline 

Typically, computations in four- dimensional gauge theories can only be carried out in 
the weak coupling regime, where a wealth of perturbative techniques have been de- 
veloped. Dualities in field theory, however, provide new avenues in which to study 
the strong coupling behavior of certain field theories by mapping the strong coupling 
dynamics of one theory to the weakly coupled regime of the dual theory. Moreover, 
some gauge theories are holographically dual to quantum gravity with certain asymp- 
totics, and the strong coupling dynamics of the gauge theory can be solved in terms of 
semiclassical gravitational physics. 

The best understood example. A/" = 4 super Yang-Mills, is a field theory that 
both presents electric-magnetic duality [H [21 E] and describes holographically quantum 
gravity with AdSij x 5*^ boundary conditions [1]. Due to its high degree of symmetry, 
it enjoys remarkable properties in the large A^ limit - such as integrability - that allow 
for the study of some questions in the strong coupling regime (see for instance [5]). 

In this paper we compute the large N, strong coupling correlation functions of a 
supersymmetric Wilson loop in a large representation of the gauge groupu with local 
operators of A/" = 4 super- Yang Mills, specifically with chiral primary operators and the 
stress tensor. We compute these correlators both in gauge theory and using the dual 
supergravity description. In gauge theory we obtain strong coupling results by solving 
the normal matrix model that captures these correlation functions. We then perform 
a quantitative test of S-duality of A/" = 4 super Yang-Mills by also calculating in the 
semiclassical approximation the correlator between a 't Hooft loop operator and the 
same local operator. We find that the S-dual of the semiclassical 't Hooft loop correla- 
tor exactly matches the strong coupling result of the Wilson loop correlator, providing 
a quantitative test of electric-magnetic duality in A/" = 4 super Yang-Mills. We also 
perform the calculation of the Wilson loop correlation functions using "bubbling" ge- 
ometries and find exact agreement with the strong coupling results we obtained in the 
gauge theory. 

The study of supersymmetric Wilson loops in the context of the AdS/CFT corre- 
spondence [H El [7] is important for several reasons. Among them is the fact that these 
operators couple to strings and branes in the bulk, thus touching on stringy properties 
of the theory. Moreover, Wilson loops allow in some cases, for instance when they 
follow circular contours, to obtain results that are exact in A^ and the 't Hooft cou- 
pling A = Qym^- The first example of an all order computation in A^ and A was the 
computation of the expectation value of the circular Wilson loop in the fundamental 



^This is a representation where the number of boxes in each row or column of the corresponding 
Young tableau is of order TV. 
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Figure 1: We depict here, rotated and inverted, the Young tableau R of the irreducible 
representation of U{N) in which we take the Wilson loop operator. The tableau consists 
of g blocks, the J-th one of them having rij rows of length Kj. All the edges of the 
diagram are taken to be long, meaning that nj and Kj are both of order 0{N) for all 
/. This guarantees that the dual bubbling geometry has small curvature everywhere. 



representation, which was conjectured in P, [9] to be captured by a hermitian matrix 
model. This result has recently been proven using localization techniques [10] and 
generalizations thereof have been found in, e.g., [HI [T21 HSl HH [15] • 

It is a well-established entry in the AdS/CFT dictionary that a supersymmetric 
Wilson loop in the fundamental representation corresponds in the bulk to a classical 
string surface with AdS2 induced metric, which extends in the interior of the AdS space 
and lands on the loop on the boundary [161 [H] • I^i particular, the expectation value 
of the loop operator is given by the minimal area of the string world-sheet, upon the 
appropriate regularization of the IR divergence associated with the infinite area of the 
string [18j . 



When one considers Wilson loop operators in representations higher than the fun- 
damental, with rank of order 0{N), the string in the bulk gets replaced by configura- 
tions of probe branes with electric flux dissolved in their world- volumes [12] (see also 
[201 |2Tl [221 [23] ) . More specifically, a Wilson loop in the rank k symmetric representation 
is described by a D3-brane with k units of flux and wrapping an AdS2 x S"^ subspace 
[201 122] , whereas a loop in the rank k antisymmetric representation is described by a 
D5-brane [ISl [21], also with k units of flux, wrapping an AdS2 x S^ subspace. These 
branes can be thought of as emerging from k coincident strings via the Myers polariza- 
tion effect [21], which, for large enough fc, blows up an S"^ or an 5^ from the world-sheet 
of the coincident strings. 

This probe approximation breaks down when the representation R of the Wilson 
loop is taken to be even larger, with a corresponding Young tableau containing order 



0{N'^) boxes (see Figured]). In this case the back-reaction of the brane configuration in 
[T9] cannot be ignored and the space-time is deformed into a new geometry containing 
bubbles of new cycles carrying fiuxes, and it is thus called a bubbling solution. The 
study of the bubbling solutions for this case was initiated in [25l EH] (see also [27]) 
and culminated in [25] , where an elegant description of the solution was given in terms 
of a Riemann surface. If the representation of the Wilson loop is large, then the 
dual geometry is guaranteed to have small curvature everywhere. The matrix model 
prediction for the on-shell action of these bubbling solutions was found in [29] . 

This entire picture of strings metamorphosing into branes and branes undergoing 
geometric transitions into new geometries is remarkably universal in the context of 
gauge theory /gravity dualities, having a very close analog for local operators, where 
these phenomena were in fact analyzed first [SUl I31]p cind for topological theories 
[331 [Ml [35]. 

The aim of this paper is to deepen our understanding of Wilson loops in large 
representations of the gauge group, both from a gauge theory perspective and in the 
bulk. The tool we use to do this is the detailed study of the correlation functions of 
these loops with local operators of A/" = 4 super Yang-Mills, such as chiral primary 
operators and the stress tensor. We perform computations both in field theory and in 
super gravity. 

In the gauge theory, the computation of these correlation functions makes use of a 
matrix model (the computation of the correlators when the Wilson loop is in the fun- 
damental, symmetric and antisymmetric representation was performed in [361 l37])n 
Building on the results in [ID], we solve the matrix model capturing these correla- 
tion functions in the strong coupling regime and when the Wilson loop is in a large 
representation. 

On the other hand, in the bulk analysis we use holographic methods on the bubbling 
supergravity backgrounds to extract the desired correlation functions, finding perfect 
agreement with our computations in gauge theory. The bubbling supergravity back- 
grounds indeed contain non-trivial dynamical information about correlation functions 
(see [¥l] |l2] for the computation of correlation functions from bubbling geometries for 
local operators and surface operators respectively)]^ 

These Wilson loops in large representations and their dual bubbling geometries rep- 
resent an arena with an incredibly rich structure, yet one where explicit computations 



^The bubbling construction for chiral primary operators was found building on ideas introduced 
in [35]. 

•^For the corresponding Wilson loop expectation value computation see [38 l l39 l [29] . 

^See [43l|37] for the probe string/brane computation of the correlation functions when the Wilson 
loop is in the fundamental, symmetric and antisymmetric representation. 



and highly non-trivial quantitative tests of the AdS/CFT correspondence are possible 
and where one may be able to shed new light into the inner workings of holography. 

1.1 Outline of the paper 

In Section [21 we establish the notation and analyze the general structure of the correla- 
tion functions we are interested in. The position dependence, both for correlators with 
chiral primary operators [3S] and correlators with the stress tensor |31|, is completely 
determined by the symmetry of the system, so that the computation boils down to 
finding coefficients which depend on the characteristic data of the Wilson loop and the 
local operator, as well as A and A^. We show using supersymmetric Ward identities 
that the correlator of the Wilson loop with the stress tensor can be obtained from 
the correlator of the Wilson loop with the dimension two chiral primary operator. In 
Appendix [B] we derive the same relation using a topological field theory argument. 

In Section [31 we proceed to compute in gauge theory the correlation coefficients 
mentioned above. The important point to stress is that we manage to perform these 
computations at strong coupling, making it possible to compare and match them with 
the results in supergravity of Section [H 

More specifically, we begin in Section 13.11 with the computation of the correlator 
between a half-BPS circular Wilson loop and a chiral primary operator of A/" = 4 
super Yang-Mills as well as with the stress tensor. Similarly to what happens for 
the expectation values, it has been conjectured in [32] that the exact path integral 
describing this correlation function is also captured by a matrix model, which sums 
all the ladder/rainbow diagrams in the perturbative expansion and therefore allows to 
extract its strong coupling behavior (a derivation using localization along the lines of 
PH] should also be possible). The particular matrix model we use is the normal matrix 
model introduced in this context in [33] ■ Using results obtained in [ID], we solve the 
model for large 't Hooft coupling and large representations of the gauge group and 
find the moments of the eigenvalue distribution, in terms of which the correlators are 
determined. 

We then calculate in Section 13.21 the strong coupling correlator between the same 
half-BPS circular Wilson loop and the stress tensor of A/" = 4 super Yang-Mills using a 
semiclassical computation of the correlator of the 't Hooft loop with the stress tensor. 
This consists in computing first the correlator of a half-BPS 't Hooft loop operator with 
the stress tensor in the semiclassical approximation. To obtain the strong coupling 
result for the Wilson loop we act with S-duality on the 't Hooft loop semiclassical 
result. This computation yields precisely the same answer computed by the matrix 
model in section 13.11 for a specific choice of representation of the Wilson loop. We 



comment on the reason why this happens. This yields a quantitative test of S-duahty 
in A/" = 4 super Yang- Mills. 

The supergravity analysis of these same correlators is contained in Section HI First, 
we briefly review the bubbling solution found for the first time in closed form in [28] 
and re-express it in terms of the resolvent of the matrix model, which encodes the 
correlation functions in the gauge theory. We then apply the Kaluza-Klein holography 
machinery |l5l HI] to this geometry and extract from the asymptotic expansion of the 
supergravity fields the correlation functions with chiral primary operators and the stress 
tensor (see also [HI US]). Differently from the gauge theory computation, we are able 
to carry on the bulk computation only for operators up to dimension four, albeit the 
procedure we use is in principle applicable to operators of arbitrarily high dimension. 
The correlators that we compute in supergravity are found in perfect agreement with 
the strong coupling gauge theory results. Particularly remarkable is the agreement 
between the correlators with dimension four operators (both the dimension four chiral 
primary operator and the stress tensor), because of very delicate cancellations between 
non-linear terms that take place in supergravity, as expected from the strong coupling 
gauge theory analysis. 

We conclude the paper with a series of appendices in which we collect some technical 
details of our calculations. 



2 Symmetry analysis of the Wilson loop correlators 

In this paper we study the half-BPS circular Wilson loop of TV = 4 super Yang-Mills 
in M^. It is given by 

WR{e,a) = ^^Tr^Pexp / {tA + <f)¥\x\ds) . (1) 

dim it J circle 

The trace is taken over an irreducible representation R of U{N). The Wilson loop 
Wr, besides containing the holonomy of the gauge field A = A^^dx^, also couples to 
the scalars 0* of the A/" = 4 multiplet through 6^, a constant unit vector on M^. The 
integral is taken along a circle of radius a in M*^ parametrized by < s < 27r. 

The circular Wilson loop is related to the straight Wilson line 

PF«"^(0) = ^_TrHPexp /" {iA + (l)'9'\x\ds) (2) 

dim R Jiine 

by a conformal transformation (an inversion around the origin). Despite this relation, 
the expectation value of the straight Wilson line is trivial (independent of the 't Hooft 



coupling A = Qyu^ ^^^ A^), while the expectation value of the circular Wilson loop 
has a non-trivial dependence on A and iV |8l |9] , which can be interpreted as a conformal 
anomaly [9]. 

We shall see in the following that it is also useful to study the physics of the half- 
BPS Wilson loop of A/" = 4 super Yang- Mills by considering the theory on AdS2 x S*^ 
instead of M^ [S]. In this geometry the Wilson loop is inserted along the boundary 
of A(iS'2lj The straight Wilson line corresponds to taking the metric on ^^52 in 
Poincare coordinates while the circular Wilson loop corresponds to taking AdS2 in 
global coordinatesO As we show in Appendix |Al the metric of AdS2 x S'^ in global and 
Poincare coordinates is related to the metric in M^ by a Weyl transformation, which 
allows us to relate the computations in M^ to the computations in AdS2 x 5*^. 

We now proceed to study the correlators of the half-BPS Wilson loop with chiral 
primary operators and with the stress tensor in A/" = 4 super Yang-Mills. 

2.1 Correlators with chiral primary operators 

A chiral primary operator of dimension J is a scalar operator transforming in the 
representation of the 5*0(6) R-symmetry with the highest weight (0, J, 0). It is given 

bjB 

where Y = Ci^,,,ij9^^ . . .9^-' is an SO{Q) scalar spherical harmonic that specifies the 
operator Oj. We normalize the complex coefficients Ci^,,,ij so that 

corresponding to the normalization of the spherical harmonics Y given by 

The operators ([3]) are normalized such that their two-point function is unit normalized 
in the planar approximation 

{Oj{x)Oj{y)) = ^ ■ (6) 

F — 1/1 



Ojix) = L_2_Q^ ,^Tr0^i . . .0^^(x) , (3) 



'''To be precise, the Wilson loop should be inserted on a curve at a finite distance away from the 
boundary and then we should take the limit in which this curve approaches the boundary. We refer 
the reader to [3? for details on the procedure. 

^Throughout this paper AdS refers always to Euclidean AdS space. 

^Here (pi = 0fT°. The gauge group generators are canonically normalized by Tr(T°T'') = i^"'' 

and the scalar propagator is given by ((/)°(a;)(/i5(0))yM = ^r "^^ . 
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(9) 



Let us first consider the correlator of tlie Wilson loop and Oj{x) in AdS2 x S"^. Be- 
cause this space is homogeneous (all points are related to each other by an isometry), 
the correlator is independent of the position x. By taking into account the transfor- 
mation properties under the R-symmetry group, we can parametrize the correlator 
as 

{Ojix))w = ^R,jY{e). (7) 

The notation {...)w denotes a correlator of A/" = 4 super Yang-Mills in AdS2 x S^ with 
the Wilson loop Wr{9) on the boundary, normalized so that {l)w = 1- This expression 
holds both when AdS2 is in global and Poincare coordinates, corresponding to inserting 
a circular Wilson loop and straight Wilson line respectively. 

For the circular Wilson loop ([1]) in M^, the correlator is given by 

{Wnie,a)Ojix)) _^ 1 

{WR{e,a)} r-' 

where we have defined the conformally invariant distance r as 

^_ ^(r2 + L2 - a2)2 + 4a2L2 

which is also the conformal factor that relates the metric on M^ to the metric on 
74^5*2 X S"^ in global coordinates, as we show in Appendix El Here a is the radius of 
the S*^, L is the distance between the location of the chiral primary operator ([3]) and 
the plane that contains the circle. The other parameter r is the distance between the 
location of the operator ([3]) and the axis of the circle. Therefore, symmetries determine 
the correlator between a chiral primary operator and a circular Wilson loop up to the 
coefficient S^j^j. 

Similarly, the correlator of the straight Wilson line in M^ ([2]) and the chiral primary 
([3]) is given by 

{W'r{e)Oj{x)) = En,jY{e)^, (10) 

where / can be interpreted, again, both as the distance between the line and Oj{x) as 
well as the conformal factor relating the metric on M^ to the metric on AdS2 x S"^ in 
Poincare coordinates (see Appendix El) . 

Agreement between various computations we perform in this paper suggests that 
the correlator of a Wilson loop with a local operator normalized by the Wilson loop 
expectation value transforms simply under conformal transformations, so that H^ j 
is the same for the correlator computed with the circular Wilson loop and with the 

8 



straight Wilson line. The transformation properties of this ratio seem to be solely 
determined by the representation of the local operator under the conformal group, 
and does not suffer from the conformal anomaly of the Wilson loop expectation value 
[9]. It would be interesting to prove this lack of conformal anomaly of the normalized 
correlator from first principles. 

The coefficient S^ j depends on R, the representation of the Wilson loop, and J, 
the dimension of the chiral primary operator (it is also a non-trivial function of the 't 
Hooft coupling A and N) but is independent of the choice of operator in the (0, J, 0) 
multiplet of 5*0(6) one uses. We will take advantage of this by choosing convenient 
operators in the multiplet for the various computations. 

2.2 Correlator with the stress tensor 

The correlator of a half-BPS Wilson loop with the stress tensor is also essentially fixed 
by symmetries up to a scaling function h\Y, that we wish to compute. For the stress 
tensor one must take into account that U{N) A/" = 4 super Yang-Mills on a curved 
background has a conformal anomaly given by 

(T„") = 3^(fl,„fl--3fl^), (11) 

where i?^jy is the Ricci tensor and R is the Ricci scalar of the background. This relation 
is protected from quantum corrections because the superconformal symmetry relates 
the stress tensor to derivatives of the R-current. 

On the AdS2 x S"^ background we are considering (both in global and Poincare 
coordinates), the anomaly flTTj) is N'^/Sn'^. This, together with the symmetries of the 
problem, determines the correlator up to a real number h^ - the scaling function - 
that depends on gyM, N, and R, but not on 9 [44J 

{T^u{x))wdx''dx'' = hw{ds\^s2 - ^40 + ^^(^4d52 + ds%). (12) 

We now turn to A/" = 4 super Yang-Mills in M^. In this case there is no conformal 
anomaly. The correlator of the straight Wilson line with the stress tensor in M^ is given 
byHH 

{W'r{e)T,,{x)) = ^ , (l^r (^)T4.(x)) = , 



{Wr{e)T^,{x)) = -h j^' ,f°^^ , (13) 



where we have taken the hne along the 4-direction and n"" = x"'/l for a, 6 = 1,2,3 is 
the unit normal vector to the line {i.e. n°'n°' = 1). 

The form of the correlator between the circular loop and the stress tensor in M^ 
can be obtained similarly and we write it for completeness. The circular loop in the 
coordinate system 

dsl^ = dr^ + Pdip^ + dL"^ + L'^dcf)'^ (14) 

is supported at r = a and L = 0. The correlator is then given by 

{Wn{e,a)Trr{x)) , /I 2r'L^\ {WR{9,a)T^4x)) , r^ 
nw \ ^ ^p^ 5 ,TT. //I — T^ iT-w^ , 



{WR{e,a)TLL{x)) _ (I {a' + L^-r^f \ {WR{e,a)T^^{x)) _ L^ 

{WR{e,a)) ^Vf4 2a2f6 J' {WR{e,a)) ^ F^ ' 

{WR{e,a)TrL{x)) , rL{a^ + L' - r^f ,_.. 

{WRie,a)) --'^ ^^ • ^''^ 

It is completely determined by the scaling function h^. 



2.3 Correlator with the stress tensor from Ward identities 

We now wish to derive an exact relation between the Wilson loop correlator with the 
stress tensor and the Wilson loop correlator with the dimension two chiral primary 
operator. This will allow us to compute the first correlator from the knowledge of the 
second one. 

The relation between these two correlators is a consequence of a Ward identity. 
The idea is to apply the supersymmetry Ward identity to operators in the current 
supermultiplet, which contains both the dimension two chiral primary operator and 
the stress tensor 



Under the Poincare supersymmetry transformations of A/" = 4 super Yang-Millq 

6A^ = -i(A^a^r7^ + r^^a^A^), (16) 

the straight Wilson line 

W'tiO) = ^T^Tr^P exp t J ds{A, + 9'E\^^^'') (17) 

is invariant under the following supersymmetries 

r]^ = te^j:\^a,v^, (18) 

^We do not write the four-dimensional spinor indices for clarity. 
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where S^^, E are the six-dimensional chiral sigma matrices satisfying (see e.g. [17]) 

S^^S^'^^ + S^j.^r^'^ = 2,5^^^, (19) 

r-^^ = -ie^^^^E^^^, (20) 

and the scalars in the 6 of SU{A) are given by 0* = S^^cy^"^^. One can regard S^^ as 
the Clebsch-Gordan coefficients couphng two 4's to a 6 of SU{A). 

Let's now consider the following correlator 

{w'rmo{x)), (21) 

where 5 denotes a supersymmetry transformation generated by the supersymnietries 
preserved by the Wilson line fITS]) and 0{x) is an arbitrary local operator. Since 
SWl'^^e) = we have that 

{W'r{0)5O{x)) = {6 {W'r{0)O{x))) = ([r/-Q,, + rJ^Ql,iy]j-(^)0(x)]) =0 (22) 

for any local operator 0{x). 

The supersymmetry variations of the supercurrent J^a (in the 4 of S'f/(4)) and the 
fermionic operator x^ab = ~x'~'ba (i^ the 20 of SU{A)) in the current supermultiplet 
are given by (see e.g. [IS] ) 

SJf,A = -(^"T^uVa - '^{(^p^^iu - -<Jf,^ap)d'"RP'^AVc 

1 

3" 



- (^pafT^. + 77^/..^P^)eACDi?5"5^^^'^r/^ 



'^X^'ab 



4 



i^ABEF(^ ijpv V + ^^ABEFC- "T] 



- a^R^yXHB] + Sza^'a^g^^A^r/B 



trace , (23) 



where a^^— trace = a^^ — (l/3)(5ja^^— 5^a|^^). The supersymmetry transformations 
generate other operators in the current supermultiplet. For example, in the right hand 
side of ([23D we get 



•\AB 1 v^^-^x^ ^„ I AiAJ xiJAkik 



QcD = 4^S, S,cDTr (^0V - -^^V^r J , (24) 

which is the dimension two chiral primary operator in the 20' of SU{A). The R- 
symmetry current R^\, the scalar operator S^^ = £^^ and the two- form -BdP = 
—B^j^ transform in the 15, 10, 6 representations of SU{4) respectively. 

We can now constrain the correlator of the straight Wilson line W^^'^iO) with these 
operators by using the fact that ^^^^{9) is 5*0(5) invariant. Since the 15 and 10 
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representations of SU{4) do not contain an 5*0 (5) singlet in the decomposition of 
S0{5) C SU{4:) we have that 

(W^|'^^(^)i?^^s) = , {W^^''%e)S^^) = . (25) 

On the other hand, since 6 -^ 1 © 5 under the decomposition, we have that 

(W^r (^)i?r) = 0, (W^r (^)5r) = b^fe^'-^ . (26) 

Likewise, we have from (ITOl) and (|24|) that 

{W'rim^cn) = ^^R,2^f^jCD{e'e^ - l^^ , (27) 

and, as we have aheady seen, 

{Wr{e)T^,ix)) = -hjj^i^p^ . (28) 



These correlators are completely characterized by the functions b, hw and S/j2, that 
depend on the representation i?, on Qym and N . 

Our goal is to relate these three quantities. For that we use the Ward identity (!22|) 
and the supersymmetry transformations (123|) . To do that, let us compute 

= {w't{e)5j^A) 



-O- -r^riA - (0-cd0-4a + -(y^a^ycdYACDEbecdeff' ( — ) 6''S T] 



3 ' I 
= (IVr (^)^x'ab) 
_ 3 

~ 4 



^W - T;b) i ^ (29) 






+ ^la-^-En^^Y^T^.ABiQ'Q' - \^'')da {^ V^ 



trace 



= l{b + ^^^R,2)^0'{^.abS^o - Is^B^^D + 1&bd)v'' , (30) 

where we have used that the supersymmetry transformation is generated by a spinor 
satisfying flTSl) . Therefore we obtain that 
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This relation has been checked at weak couphng to make sure that the numerical 
coefficient is correct. We stress that the relation holds exactly for arbitrary R, gyMi 
and N , as it follows from a Ward identity. This allows us to calculate the correlator 
of the half-BPS Wilson loop with the stress tensor in terms of the correlator of the 
Wilson loop with the dimension two chiral primary operator]^ This will allow us to 
compute the stress tensor correlator at strong coupling by solving a matrix model. 

We expect that similar arguments can be constructed to relate the correlator with 
stress tensor to the correlator with the dimension two chiral primary, in cases involving 
other supersymmetric operators/backgrounds, e.g. surface operators, half-BPS local 
operators, and interface CFT's. 



3 Correlators from gauge theory 

In this section we compute the coefficients S^ j and hw in field theory. Later, in Section 
m we will repeat these computations in supergravity using the bubbling Wilson loop 
supergravity solutions. As explained earlier, 'Er^j is independent of the choice of Cj^..^^, 
i.e., the choice of the spherical harmonic Y{9)., or, equivalently, it is independent of 
the choice of operator in the (0, J, 0) S'f/(4) multiplet. We will take advantage of this 
when calculating S^j j in gauge theory and supergravity. 

3.1 Correlators from a matrix model 

3.1.1 Complex and normal matrix models 

So far we have not committed to any explicit choice of chiral primary operator repre- 
sentative in (l3l). We do this now in order to compute S/j^j, having in mind that the 
final result is in fact independent of this choice. We take the following definition of 
complex scalar field: 

Z.^, (32) 

and choose the following chiral primary operator: 

Oj(x)= ^ ' TiZ-^, (33) 



®See Appendix |B] for an alternative derivation of this relation obtained using a topological field 
theory argument based on the GL twist ^49J . 
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which, as we shall see shortly, allows for its correlator to be computed using a matrix 
model. The 9 dependence of the correlator of this operator with the circular Wilson 
loop Wr{6, a) is given by ([8]) with 

yw - <q#*^. (34) 

Using the symmetries of the problem, we can specialize without loss of generality 
to a configuration with r = 0, corresponding to the local operator being inserted on 
the symmetry axis of the circle. We can moreover use the 5*0(6) symmetry to take 
6 = (1,0, ...,0) so that the Wilson loop ([T]) only couples to 0^. In this case the 
contribution to the correlator ([8]) due to the spherical harmonic associated with (1331) 

is Y{e) = 2-^/2. 

It was conjectured in ^36j that radiative corrections to this correlator that involve 
internal vertices cancel to all orders in perturbation theory and therefore do not con- 
tribute to the evaluation of correlators between chiral primaries and circular Wilson 
loops. This is the working assumption we make for the gauge theory analysis (a deriva- 
tion using localization similar to the one in [TU] should be possible). Moreover, with 
the choice r = 0, every point on the circle is equidistant from x and the propagator 
between the chiral primary and the Wilson loop becomes constant. It was first noticed 
in [8] that, in Feynman gauge, the combined propagator for the gauge field and the 
scalars between two points on the circle is also a constant (independent of the radius 
a of the circle). 

Summing over all Feynman diagrams reduces then to a combinatorial problem, 
where one has to count the number of free propagators at any order in perturbation 
theory. As pointed out in [3H], this combinatorics is exactly captured by a complex 
Gaussian matrix model defined by the partition function Zc = f[dz] exp (— ^TrJz), 
where z is a complex N x N matrix. This matrix model also computes the two-point 
function of local operators in A/" = 4 super Yang-Mills of the form Tt{Z'^) [50]. There- 
fore, the correlator of the circular Wilson loop ([T]) with the chiral primary operator 
( l33l) is given by 

{WR{e,a))YM rJ XJ/^y/j /[rf2]e-^Tr^^Trfie(-+^)/2 ' 

where r is given in iQ . 

By comparing this expression with ([8]), and using that for 6 = (1,0,..., 0) the 
spherical harmonic function corresponding to fl33l) is given by Y{9) = 2~'^/2, we find 
that 

2^/2 J[dz]e-'-^^'^'TiRe('+^y^Tiz-^ 



-Rj 
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(36) 



Therefore, we arrive at the result that the correlator of a dimension J chiral primary 
operator Oj with a half-BPS circular Wilson loop Wji{d, a) is captured by the moment 
{z'^) of a Gaussian complex matrix model fl35|) . 

After having carefully settled the normalization factors, the next task is to compute 
the moments {z'^) in the complex matrix model. This is not easy because the eigen- 
values of the complex matrix model do not decouple from the off-diagonal components 
of the matrix. To proceed, we will map the complex matrix model to a normal matrix 
model, where one can reduce the computation of the moments to integrals over the 
eigenvalues. Via manipulations involving coherent states, the authors of [38] proved 
the following formula (see Appendix A of [38j): 



1 



i 

= ±. I <f^ze-'^^^-^^^^A{z)\'\{e'^'^e-^''l (37) 

On the left hand side we have the hermitian matrix model, with eigenvalues ^j and parti- 
tion function Zh = J d^^e-^^^^^ /\{if, while Zi and Zm = / rf^iv^g-Sf e,^.-.|A(^)|2 
are the eigenvalues and partition function of the normal matrix model. The factors of 
A are the Vandermonde determinants originating from the transformation to the eigen- 
value basis and the constants ki encode all the information about the representation 
of the Wilson loop insertion. 

Since Trij(e^) is a polynomial of e^% the equation above proves that the her- 
mitian and the normal matrix model are almost equivalent, upon the replacement 
Tr/j(e^) -^ Tr/j(e''^"*"^^/^^). Because of the extra factors e~»N^t in the right hand side, 
this equivalence seems limited to the anti-symmetric representations where these fac- 
tors are independent of the index i and can be pulled out of the integral [38] . We can 
circumvent this difficulty by rewriting 



e-w^? = V ^ / ^«^e-^"H^fc.".. (38) 

Thus we find that (going back to the matrix form for conciseness) 

= ^l Mzl[<ia]e-^-^P^+'>")-i-Tr„(e«+">), (39) 

^N^a J[z,z]=o dim it 

where [d^] and [dz\ are the hermitian and normal matrix measures, respectively, a is 
a real diagonal matrix, and [da] is the Euclidean measure. We have divided by the 
appropriate partition functions for proper normalization. 
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We are in fact interested in an extended version of the relation above, which is 
obtained by applying the trick fl38|l to the results discussed in Appendix C of [38] . 
This extension includes an extra insertion corresponding to a chiral primary operator, 
and is given by 



= -L-i- / HMa]e-^-(--"^)-i-Tr^(e^--)Tr.^ 
2JI^ Zf^Z^ Jy^^^^Q dimi? 

^ l,WnTxz')MM (40) 



2J/2 



This formula rewrites the complex matrix model correlator in (!36|) as a normal ma- 
trix model correlator. We will calculate H/j^j in (136|) using the normal matrix model 
description. 

We are interested in the computation when the representation R is large while J 
is of order one. It is for this class of operators that the dual supergravity background 
that we will compute with in the next section has small curvature everywhere. In this 
case, the eigenvalue distribution of the hermitian and normal matrix models is altered 
by the Wilson loop insertion but not by the chiral primary operator insertion. 

The trick found in [35] and used in [40j to analyze the eigenvalue distribution of a 
hermitian matrix model with a Wilson loop insertion in a large representation R (see 
Figure [1]) is to split the matrix ^ into (7 + 1 blocks ^i of size ni x ni and rewrite the 
traces in terms of interactions among different submatrices 



dim/?(W^^)MM = ^ [me-'-^^^'TiRie- 



2 



= ^ / wu e-^^.-^h^.-'-'-^' n d<'t '^; ^ '~'^ f:' . (41) 

Zh J ^f^ f-^^, 1 - e ^10 e^i' 

where Kj are defined in Figure [1] At the saddle point of the integral, the eigenvalues 
of C,i for fixed / are distributed along some interval [e2/,e2/_i]. These intervals are 
ordered as 

629+2 < . . . < ei . (42) 

In the limit 

A»l, g'y^jnj = OiX) , gl^,{Ki - Ki^,) = 0{\^/^) , (43) 

these cuts are separated from each other by a distance of order vA and the exponential 
interactions in (HTl) can be safely ignored. We note that we are studying the matrix 
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model in the supergravity regime fH3l) . where it is meaningful to compare the matrix 
model computation with the corresponding computation performed using the dual 
bubbling supergravity solutions, which we carry out in the next section. Thus the 
matrix integral for the Wilson loop expectation value in a large representation R (see 
Figure [T]) in the supergravity regime fH51) is given by 

(XuivR{Wr)mm 

=^ I w^^^ ''"^^^ ''''''^' '''^''' n det (0 ® 1 - 1 ® ijf ■ (44) 

^ -^ I KJ 

This model was solved in [30] at large A^. 

The same trick goes through also in the normal matrix model. It is also straight- 
forward to generalize the computation with a local operator insertion, which is what 
we need to compute (HUl) . The corresponding formula is given by 



dim R{WrTtz'^) 



MM 



J[zi,zi]=0 J 

n Idetfe«l-lg'^^-)l' ^ V TV^/. (45) 



ZfqZa J^zi^i]=f) 



X 

i<i' det ( 1 - e"^7r-'"-f (g, e^TJ 






Similarly to the hermitian matrix model, in the supergravity regime (H3|) . the expo- 
nential interactions in the denominator in the second line of ( l45ll can be neglected. 
Then the dynamics of a is decoupled and trivial as it doesn't have the Vandermonde 
repulsion. Thus we get 



{WrYiz') 



MM 



^b)mm 
f ll[dz,]e~"-^^'^^^'^'e^''''^('^) H \det{zi^l-l^zi,)\'J2^Tz/ 

I \{[dzi]e~"-^^^"'^-^''''e^^'''^i'^) W \det{zi ® 1 - 1 ® z,,)!' 

(46) 

We can study the resulting matrix model in the large N saddle point approxima- 
tion. Since the submatrices of z feel different constant forces proportional to Kj, the 
eigenvalues of z spread into g + I droplets (see Figure [2](c)), and we can easily obtain 
the normal matrix model saddle point equations 

2N Ki ^^ 1 
---^7i + -=+ > = ^ = 0, 1 = 1,. ..g + 1, z = l,...,n,, (47) 
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where / labels the droplets Dj and i the eigenvalues inside each droplet, and Yl^i=i ^i = 
A^ (see Figure [1]). 

3.1.2 Large A^ solution of the normal matrix model 

In order to solve the normal matrix model in the supergravity regime fH3|) . we define 
the resolvent of the normal matrix model as 



T^ z- zj^i Jc z- z' 



(4J 



where (t{z^^) is the eigenvalue density in the complex plane and we use the measure 
(Pz = d{R.ez)d{lYaz). As in any normal matrix model, the eigenvalue density is con- 
stant in the droplets (which are then incompressible), as one easily sees by rewriting 
the large A^ saddle point equations (H7|) as 

-2z+^^^ + n(z) = (49) 

v2 

for z G Dj, and using that dz[l/{z — z')] = nd^i^z — z'): 

^ for z G Dj, 

(50) 
for z ^ Dj. 



The resolvent then becomes 



n(z) = - f d^z' —^ . (51) 

t^Jd=\JjDj z-z' 

Our task is to find a function VL{z) and g + 1 simply connected regions Dj such that 
VL{z) is holomorphic outside D = UiDj, equation fj^9|) is satisfied on the boundaries 
dDj, and 

n{z) = - + 0{z'^) asz-^oo. (52) 

z 



As we saw in equation ( 139I1 . the normal matrix model is simply a rewriting of the 
hermitian matrix model if we ignore the Trz'^ insertion, to which the eigenvalues do 
not back-react anyway. So we should expect that the large A^ solutions of the hermitian 
and the normal matrix models are related. 

The hermitian matrix model (IHl) has been solved in HOl. The resolvent 



^1 (C) = 9ym Y1 737C0 ^^^^ 
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(a) 



(b) 






Figure 2: (a) The branch cuts [e2/, e2/_i] (/ = 1, . . . ,g + I) on the ( plane C for the 
the hypereUiptic curve ( l56l) in the g = 2 case, (b) The eigenvalue density p(^) of the 
hermitian matrix model ( 1411) . (c) The corresponding eigenvalue distribution (droplets) 
of the normal matrix model (l45l) . The shape of a droplet is given by p{C,)- 



is given as the indefinite integral 



^i(C) 



c 



s+i(C) 



Here ag+i{() and 



2-2- 

oo \ ^JH2g+2{C') 



2g+2 
H2g+2{()= n(C-e. 



rfC'. 



(54) 



(55) 



are monic polynomials of degree g + I and 2g + 2, respectively. Their coefficients are 
determined by the constraints described in [lU]- These constraints guarantee that the 
integration contour in fl5^ is arbitrary as long as it does not cross any of the cuts 
[g2I, e2/-i]. The integrand is a meromorphic one-form on the hypereUiptic curve given 
by the equation 



w 



H2g+2{0- 



(56) 



Let us define a holomorphic function uj2{C) by analytically continuing the resolvent 
co'i(C) to the second sheet along the {g + l)-th branch cut: 



U2{^ ± ie) = u;i(^ =F ie) for ^ G [62^+2, 62^+1] 



In other words. 



C 



«3+l(C') 



^2(0 = ^1(629+2) + / 12+ 

'e2g+2 \ V-"23+2l<; ) 



dC, 



(57) 



(5J 
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where again the contour should not cross any cut. It then follows that 

:,, (59) 






where the contour Bj goes around the interval [623+1,62/] , from 62/ to 62^+1 on the 
first sheet and from 62^+1 to 62/ on the second sheet. The last equality in fl59|) is one 
of the constraints HO] . 



We claim that the solution to our complex analysis problem determining the resol- 
vent of the normal matrix model is given by ri{z) = (l/-\/2)^i o uj2^{2^/2z), or 

^2(0, (60) 



2V2 



n{z) = ^^i(C). (61) 

First, one can see from flSSl) that flBU]) maps C — U/[62/, 627-1] to the complement of 
some droplets, identified with Dj. That is, it maps Figure [2|^a) to Figure [2](c). For 
z G dDj with Im2; ^ 0, then 

-2. + ^l:^ + ^]M=-i=c.2(e±^e) + ^i^ + i= (62) 



where we used that uJi{^ ± ie) = uji{^^ it). Thus ( H9!) is indeed satisfied! 

In summary, we have found a normal matrix model eigenvalue distribution in terms 
of the hermitian matrix model. This now allows to calculate "^r^j in (l36l) in terms of the 
moments of the hermitian matrix model eigenvalue distribution. We will then compare 
this with the supergravity computation in the next section, where exact agreement will 
be found. 

3.1.3 Correlators as moments in the normal matrix model 

All the information about the representation R of the Wilson loop is encoded in the 
moments of the matrix model eigenvalue distribution, which for the hermitian matrix 
model are defined by 

{n^Pn^ [d^piOC (63) 



Here p(^) is the hermitian matrix model eigenvalue distribution. We want to express 
the normal matrix model moments 

(.^)^ /.-.(...)- (64) 
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in terms of the (^")'s. 

The hermitian matrix model resolvent uJi{() can be expanded in moments of the 
eigenvalue distribution 

^i(C) = A/dep(OT^ = Af;|^, (65) 

with (^) = Oo This constraint can be interpreted as arising from considering A/" = 4 
super Yang-Mills with SU{N) gauge group. 

From the SU{N) saddle point equation on the [g + l)-th cut, we have for any ( 
that |40] 

iUiiO+MO = < + 9ym\R\/N, (66) 

where \R\ is the total number of boxes in the Young tableau corresponding to the 
representation R (see Figured]). After shifting z by v^(7yj^^|_R|/4A^ in fl60l) . we get that 

z = V2C ~ ^uj^iO . (67) 

This relation can be inverted recursively using (1^ 1^71) to obtain 

z X 4A(0-A2 X{e) 

C = ^ + —^ + r- o + 



V2 2V2z AV2z^ z 



4 



v^A^ - 8v^(e')A2 + 8x/2(e')A 

.5 



^- 8.^ "^ ^0{z-\ (6J 



We can now write the resolvent of the normal matrix model in terms of the hermitian 
matrix model moments of the eigenvalue distribution by combining fl6H EHD 

j=o 

1\3 



A , 2A(a^JA2 2v^(a , 4Ar)-4A2(a + |A 

z^ z'^ 

4v^A(e')-5x/2A2(e=^) 



z z^ z^ z^ 



z^ 



8{e)X - 12(e")A2 - 6(^2)2^2 _^ ^(^2^^3 _ M4 



Z^ 



+ O {z-^) . (69) 



This allows us to express the normal matrix model moments {z-^) in terms of the 
moments of the hermitian matrix model (^") = p„. Notice that non-linearities in the 



^^ This imposes a center of mass condition on the distribution. So far we have focused on the 
U{N) case. For SU{N) gauge group, we need to restrict the integrated hermitian and complex 
matrices. Since the t/(l) part is decoupled, in the large N limit, the net effect is to shift the eigenvalue 
distributions so that the average eigenvalues vanish: (^) = (z) = 0. 
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moments start to appear only at order 0{z ^). The first few moments of tlie normal 
matrix model are given in terms of the moments of the hermitian matrix model by 

{z') = 2{e) - ^ = 2Ap, , 

{z') = 2V2{e) = 2V2Aps, 

{z') = 4^) - A\{e) + Y = 4 (Ap4 - AAp2) , (70) 

where Ap„ = pn — Pn- Here we have used that the moments p^ of the Wigner semi- 
circle distribution law, which has the eigenvalue distribution p°(0 = (2/7rA)A/A — ^^, 
are given by P2 = A/4, P3 = and p^ = A^/8. As we shall see in the next section, the 
fact that the correlators are given by moments relative to the Wigner semi-circle distri- 
bution has a corresponding statement in supergravity, where correlators are captured 
by devitations away from the AdS^ x S^ vacuum solution, which in the parametrization 
of [28] corresponds to the Wigner semi-circle distribution law. 

We are now ready to write the result of our computation of the correlator coefficients 
S/j J in the supergravity regime fll5]) . From fl5B]) . we find that 

^R,2 = V 2 — Ap2 , 

3A37^' 

N 
Sfi,4 = 2-(Ap4-AAp2). (71) 

The expression ( !7T|) . together with (jTj), is the final result of the gauge theory com- 
putation for the correlators between a half-BPS Wilson loop and the chiral primary 
operators in A/" = 4 SYM and represents a prediction for supergravity, as we have 
solved the matrix model in the supergravity regime (H5]) . 

Moreover, using our derivation of the relation between the correlator of the Wilson 
loop with the stress tensor and the correlator with the dimension two chiral primary 
operator fl31j) we also obtain the correlator with the stress tensor in terms of the 
hermitean matrix model. Namely, we have that 

A^ N^ 

and the correlator is given by flT2l) . 

In Section m we will calculate the same correlation functions using the dual bubbling 
supergravity solutions and show that we get precise agreement. 
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^i?„3 = 2\ -—j:tAps, 



3.2 Correlator of Wilson loop with the stress tensor from S- 
duality 

We have already shown in subsection 2.3 an exact relation between the correlator of the 
Wilson loop with the stress tensor and the correlator of the Wilson loop with dimension 
two chiral primary operators. Using this result, we have computed the strong coupling 
result of this correlator in terms of the hermitian matrix model second moment of the 
eigenvalue distribution (1721) . 

Here we calculate the stress tensor correlator at strong coupling in an alternative 
way. This involves considering the correlator of the stress tensor with a 't Hooft line 
jThne ^^ ^YiQ semiclassical gauge theory and then S-dualizing. This turns the 't Hooft line 
into a Wilson line and exchanges the weak coupling regime with the strong coupling 
regime. So we need to calculate the semiclassical scaling weight hr for the 't Hooft line 
T]^^*^, which captures the correlator of a 't Hooft line operator with the stress tensor 
as in (I13p (see also [B]). We find exact agreement, providing a non-trivial test of 
S-duality in A/" = 4 super Yang-Mills. 

The bosonic action of A/" = 4 super Yang-Mills is given by 

^ = ^jd'x^Tr Qf;..F^^ + D,<p'D^<p' + ^<PY - ^[0\ <P'?) ■ (73) 

The bosonic contribution to the stress tensor of the theory is given by 
2 6S 



T, 



flU 



V9^9 



liv 



2 „ /„ „ „ 1 „ \ 2 „/„,..„ ,.. 1 



Tr F,pF/ - -g,,F,,FP'^ + ^Tr D^0'D,0' - -g.^DpCJ^'D 



K" nPA,^ 



iJYM \ ^ / HYM V ^ 



(74) 



We want to compute the correlator (T]^'^'^T^,^(a;)). In the semiclassical approxima- 
tion, this is found by evaluating the stress tensor fl7^ on the gauge field configuration 
produced by the insertion of Tjj"*^, which is given by (see e.g. \M\) 

F = ^Bvols2, 0=|. (75) 

Here V0I52 is the volume form on the S'^ surrounding the 't Hooft line, / is the distance 
from the 't Hooft line and B is the highest weight vector for the representation R (see 
Figure [T]) . 
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The semiclassical scaling weight for the 't Hooft hne operator T^^'^ is given by (c.f. 

(USD) 



1 ~^^2 



3g 



2 Tt{B ) + corrections , (76) 



YM 



where the corrections are due to loop effects, that would be interesting to compute. 
For gauge group SU{N), B - the highest weight vector of the representation R - is 
given by 

B = diag (i?i - \R\/N, R2 - \R\/N, ...,Rn- \R\/N) , (77) 

so that (see Figure 1) 

TrB^ = ^ni{Ki - \R\/Nf. (78) 

Therefore, we have found that in the semiclassical approximation the scaling weight of 
the 't Hooft line operator T^^'^ in a representation R is given by 



hi 



3g 



^—J2MKi-\R\/Nr. (79) 

YM T 



We are interested in computing the scaling weight of the corresponding Wilson line 
W^^^ hw at strong coupling. S-duality is expected to exchange the 't Hooft line T^^"^ 
with the Wilson line W^]^'^'^, as well as exchange the weak coupling regime with the strong 
coupling regime. Therefore, for the Wilson loop scaling weight hy/ at strong coupling, 
we should S-dualize the 't Hooft loop result (17^ by replacing gy^^ -^ IGn^/gyM- The 
S-dual scaling weight of the 't Hooft line - which we denote by /if - is then given by 

/i| = -^^y"ni{Kj - \R\/Nf + corrections . (80) 

On the other hand, the strong coupling result we obtained using the normal matrix 
model for the scaling weight of the Wilson loop W}^^*^ is ( !72|) 

hw = -7r^^P2. (81) 



We now note that 



{z') = 2Ap, = ^E^^(^^ - 1^1/^)' (82) 
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in the limit that the cuts in the eigenvalue distribution of the matrix model are widely 
separated, so that 



,(,)^^|,(,- ^^m(a-,-W/v) ^^ 



a) 



In this particular limit, the expression for the S-dual of the 't Hooft loop scaling weight 
( IHOjl agrees precisely with the computation of the scaling weight of the Wilson loop in 
the strong coupling regime, obtained by combining (IHTl) and (!82l) . This is a non-trivial 
quantitative test of S-duality for A/" = 4 SU{N) super Yang-Mills. 

It is desirable to understand why quantum corrections are suppressed in this limit 
and to explicitly compute them. Once these are included in the 't Hooft loop compu- 
tation, the agreement we found could be extended. We hope to come back to these 
issues in the near future. 



4 Correlators from supergravity 

In this section we compute the correlation functions of a half-BPS circular Wilson loop 
with the chiral primary operator Oj for J = 2, 3, 4 and with the stress tensor of jV = 4 
super Yang- Mills using the bubbling supergravity backgrounds found in [25l 1261 l2H] ■ 
These geometries are regular solutions of ten-dimensional type IIB supergravity that 
are asymptotically AdS^ x S^ and provide the gravitational description of all half- 
BPS circular Wilson loops in A/" = 4 super Yang-Mills. They capture the complete 
backreaction of the configuration of D5 or D3 branes in AdS^ x S^ describing a half- 
BPS Wilson loop in an arbitrary representation R of the gauge group p^ (see also 

[201 EH [23). 

For such ten-dimensional asymptotically AdS^ x S^ solutions there is a well-defined 
procedure, developed in [45l|41]j^to extract the one-point functions of local operators 
in the state produced by the operator that the bubbling solution describes. Using this 
method, we will be able to obtain the correlators of the Wilson loop from the asymptotic 
expansion of various bulk fields, which we compute using the bubbling supergravity 
solution. We will find exact agreement with the strong coupling computation in gauge 
theory performed in Sections [2] and [3] using matrix models and S-duality. 



^^For some of the previous work on holographic renormahzation see [51] [52l [53l [SH [55]. A nice 
review of these topics can be found in [56j . 
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4.1 Review of the bubbling solution 

We start by briefly reviewing the bubbling solution found in [25], EH [28], using the 
elegant parametrization of the solution found in [28] (to which we refer the reader 
for more details). As is well-known, a half-BPS circular Wilson loop preserves an 
Osp(4*|4) subalgebra of the PSU{2, 2|4) algebra of symmetries of A/" = 4 super Yang- 
Mills. The 5*0(2,1) x SO{'i) x 50(5) bosonic symmetries in Osp(4*|4) are realized 
in the ten-dimensional supergravity solution by writing the ten-dimensional metric as 
an AdS2 x S'^ x S^ fibration over a two-dimensional base manifold and by writing the 
most general ansatz for the other supergravity fields compatible with this symmetry. 

The metric describing a half-BPS Wilson loop is then given by 

ds' = flds\^s, + fldsl, + flds% + Ap\dx^ + dy^) (84) 

where the warp factors /i, /2, Ja-, and p are real functions on the baseo The warp 
factors and the supergravity fluxes can be completely expressed in terms of two har- 
monic functions hi and /i2 on the base [28]. In [ID], a precise relation has been found 
between these harmonic functions and the data that control the spectral curve of the 
hermitian matrix model, which captures the vacuum expectation value of a half-BPS 
Wilson loop in A/" = 4 super Yang-Mills [H [9|, [10] . The mapping is given by [IQ] 



ia' ,^, , , ,, , ia' 



where 



hi = —{2{z-^)-{uJi-uJi)), h2 = —{z-z) (85) 

o9s 4 



ivXsmh.{x + iy) , (86) 



and uji is the hermitian matrix model resolvent introduced in equation (!53|) . and z is 
the spectral parameter of the resolvent. Notice that we use from now on, for graphical 
clarity, the letters z, x, and y for the hermitian matrix model variables (whereas in the 
previous section we have used Greek letters for the hermitian matrix model and Latin 
letters for the normal matrix model). The hermitian matrix model resolvent is given 
by 



.2 rr.^ 1 _ X /■ P(^) 



uji = g'jTi — = A / ^-^^. (87) 

yyM ^_jvf J z-x ^ ^ 

We note that the information of the representation of the Wilson loop is encoded 
in the resolvent ooi, which depends non-trivially on the choice of representation R of 



^^ The same symbol p denotes both the eigenvalue density as well as a component of the metric. 
The distinction should be clear from the context. 
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the Wilson loop, while the harmonic function /i2 is universal and independent of the 
representation. 

It is convenient to define the following combinations of the harmonic functions hi 
and /i2 and their derivatives 

_ 1 dhi dh2 1 dhi dh2 

2 dy dx 2 dx dy 

1 dhi (9/i2 1 dhi dho 

W = 1 

2 dx dx 2 dy dy 

1, , I fdhiV fdhi^^' 



hlw, 



N2 = -hih2 



dh 
dx 



+ 



dh2 
dy 



hlW. 



(88) 



The warp factors in the metric ( !84l) are then given by |28] 



f? 



fl 



L 



-^/-t4 



1/2 



-^ht^ 

N2 ^N2 



1/2 



' N2W 

W^NiN2 
h\hi 



1/2 



1/4 



^9) 



The RR four-form can be read off from 



-^4567 



dC(i) = -d]i^''' + d]2e' 



(90) 



where j2, as shown in [40], is given by 
d32 = 



■iftpifzdw - hdw) 

dx \ W J \ dy dy 

+ ^ (hih2-\ -?,(hi^ - h2 — 

dy \ W J \ dx dx 



dx 



dy. 



(91) 



In the expression for the four- form gauge field fl^Dl) . e^^^s jg ^^iq AdS2 x S'^ volume form 
and e^^^''' is the S'^ volume form, both with unit radius. 
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The complete bubbling supergravity solution also excites the dilaton and the RR 
and NS-NS two-form gauge fields of type IIB supergravity. However, we will not need 
their explicit expressions in this paper, as we will later show that these fiuxes do not 
contribute to the correlation functions we compute. Their explicit expressions can be 
found in [28j . 

In order to calculate the correlation function of a half-BPS Wilson loop and a local 
operator, we must study the deviations of the bubbling supergravity solution from the 
AdS^ X S^ vacuum solution. Our task will then be to extract the various correlation 
functions from the deviations from this vacuum. Therefore, we first consider the eigen- 
value distribution corresponding to AdS^ x S^. In the language of the hermitian matrix 
model, the eigenvalue distribution for this case is the Wigner semi-circle law 

wf^ = 2z-2Vi^^A. (92) 

From equations (l85l 1861) we get that the harmonic functions corresponding to AdS^ x 5*^ 
are given by 

/i-[ = 4c^ cosh X cos y , /i2 — 4c^sinhxcosy, (93) 

where c^ = vAa'/S (we have set the background dilaton to zero for simplicity). For 
AdS^ X S^ the functions defined in fl88|) are given by 



V^(o) = _4c4 sin 2y , 
H/(o) = 4cSinh2x, 
N^ = — 64c^ cos^ y sinh 2x , 
Arf ) = 64c^ cos^ y sinh 2x , (94) 

while the AdS^ x S^ warp factors ara^i 

/f ^ = 2V2c cosh X , fP = 2V2c sinh x , /i°^ = 2V2c cos y , p(°) = v^c . 

(95) 

These warp factors give rise to the AdS^ x 5*^ metric 

ds"^ = 8c^ (cosh^ X ds\^g^ + dx"^ + sinh^ x dsl.2) + 8c^ [dy"^ + cos^ y dsl.4) , (96) 

where AdS^ is foliated by AdS2 x 5*^ slices while the S^ is foliated by S"^'s, a slicing that 
makes manifest the symmetries of the half-BPS Wilson loop. For the vacuum solution 
only the RR four-form is excited and 

dj(°) = -64c^cos^l/t/y. (97) 



i^Herej/e [-tt/2,t:/2]. 
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4.2 Kaluza-Klein holography 

We recall now how to holographically compute the one-point functions of local oper- 
ators from asymptotically AdS^ x 5*^ supergravity geometrieso For more details we 
recommend HSl 



Given an asymptotically AdS^ x S^ supergravity solution, one needs to expand 
all the ten-dimensional fields excited in the solution in a complete basis of spherical 
harmonics on the 5*^. This produces in general an infinite number of five- dimensional 
fluctuation fields. These fluctuation modes are, however, not independent. Some of 
them are in fact related to each other by the action of ten-dimensional diffeomor- 
phisms, which give rise to non-linear gauge transformations on the five-dimensional 
fluctuations. Instead of gauge fixing these symmetries (by going for example to de 
Bonder gauge) as it is usually done in the study of the spectrum, it is more convenient 
to construct gauge invariant combinations of the fluctuations. This is because generic 
ten-dimensional supergravity solutions, such as the bubbling Wilson loop backgrounds 
we are considering, are generally not in de Bonder gauge. The equations of motion 
solved by the gauge invariant fluctuations are nevertheless the same as those of the 
fluctuations in the de Bonder gauge. 

The gauge invariant combinations of fluctuations obey in general non-linear equa- 
tions of motion containing higher derivative terms, just like the fluctuations in the de 
Bonder gauge do. These equations of motion with the higher derivatives, however, can- 
not be obtained from a local five-dimensional action. In order to perform holographic 
computations of correlators using supergravity, we want to rewrite the bulk action 
using a local bulk action. This can be accomplished by performing a Kaluza-Klein 
reduction map, a non-linear map between solutions to the ten-dimensional equations of 
motion, ipiod, and solutions to the five- dimensional ones, '^5d, which can be schemati- 
cally expressed as 

^5d = ipiod + IC ipiod ipwd + ■■ ■ (98) 

where /C is some operator containing also derivatives and the ellipses denote higher 
order combinations of the 10- dimensional fields and their derivatives. Notice that, in 
principle, all Kaluza-Klein modes are kept in the reduction map. However, in practice, 
when computing the expectation value of some operator with a given dimension, only 
a finite number of modes will contribute, giving an effective truncation of the Kaluza- 
Klein tower. This will also limit the number of non-linear terms in f l98|) that one needs 
to compute. In general, the higher the dimension of the operator one considers, the 
more terms have to be turned on. For example, we shall see that for the dimension 



^"^The procedure is actually more general and applies to any asymptotically AdSp x X"? solution, 
with X* being a compact manifold. 
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four chiral primary operators only two terms in the map are needed: a linear term with 
dimension four and a term quadratic in the fields dual to dimension two operators. 

At this point, the new equations in the reduced fields "^^d can be integrated into a 
local five-dimensional action and one can use the general holographic rules to compute 
gauge theory correlation functions from supergravity [HEllT]. The local five-dimensonal 
action does suffer from infrared divergences and has to be regularized by the addition 
of appropriate boundary counterterms and then evaluated on-shell. Differentiating the 
regularized on-shell action with respect to appropriate sources in the spirit of [U [6l [7] 
yields the renormalized one-point functions of the dual gauge theory operators. These 
are related to certain coefficients in the near boundary expansion of the bulk fields, 
which correspond to the normalizable fluctuation mode of the field. The higher the 
dimension of the operator, the deeper in AdS space one needs to dig to extract its 
one-point function. 

In this paper we are interested in computing the one-point functions of Oj for 
J = 2,3,4 and the stress tensor in the state created by a half-BPS Wilson loop. The 
calculation of these relatively low dimension operators enjoys a great simplification. In 
order to calculate their one-point functions we can neglect all the supergravity fields 
excited by the bubbling solution except for the metric and the RR four-form gauge 
field. The fluctuations that arise from the other fields (the dilaton and the RR and 
NS-NS two-form gauge fields) do not enter in the calculation of the one-point function 
of the local operators under study. Technically, the reason this occurs is that the 
fluctuations coming from the dilaton, the RR and NS-NS two-form gauge field fall off 
too fast near the AdS^ x S^ boundary and therefore do not enter into the Kaluza-Klein 
reduction map (1^ for the dual fluctuations. We can therefore use the formulas relating 
five-dimensional fluctuations and one-point functions of dual gauge theory operators 
derived in ESl ST] . 



Without further ado, we setup our computation by expanding the bubbling Wil- 
son loop solution in fluctuations around the AdS^ x S^ background. We define the 
fluctuations in the metric fl84l) as 

ds' = {fiy{l + Ai) dslds, + (/f )^(1 + A,) ds% 

+ (/f )2(1 + A4) ds% + 4(pW)2(l + A,) {dx^ + dy^). (99) 

Physical quantities, such as the correlation functions we are after, are encoded in the 
asymptotic expansion of these functions for large x. 

Since all the information of the matrix model is given by the resolvent, the fluctu- 
ations Ai, A2, A4, and Ap should depend only on the matrix model resolvent ui. So 
we first need to determine the relation between the resolvent and the asymptotic form 
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of the A's. We start by expanding hi for large x 

/ oo 

hi = 4c^ cosh X cosy ( 1 + J]] c„(i/)e-™^ ) , (100) 



m=3 



whereas h2, being independent of cui, remains equal to its background value /i2 corre- 
sponding to AdS^ X 5*^. The corrections start from order e~^^ because the background 
(coshx) includes e^' and e~^ and the corrections should not affect these terms, as the 
solution should match on to AdS^ x S^ asymptotically. The first coefficient, 03(1/), does 
not actually have any physical information because it can be eliminated by going to 
the "center of mass" coordinate system or, in other words, we can always set the first 
moment of the matrix model eigenvalue distribution to zero without loss of generality. 
As we noted in footnote [101 this is automatic for SU{N) gauge group. 

The Laurent expansion of ui in z is given by 



-. = ^Ep5t. (101) 

where pn are the moments of the eigenvalue distribution introduced in f l63|) : 

rfxp(x)x". (102) 

We would like to express the functions that appear in the asymptotic expansion of the 
harmonic function hi (IIOOI) in terms of the PnS, which contain the information about 
the eigenvalue distribution of the matrix model. This can be systematically computed 
by plugging (llOip into the definition of hi in (185|) and comparing the result with (IIOOI) . 
We get for the first few coefficients 

C4{y) = ^-^(1-2 cos 2?/), 

64Ap3 . . 

^5^^) = ^^575- COS 2?/ sm?/, 

ced/) = ^/^ (l-2cos2|/ + 2cos%) -^(2 - 4cos2y + 3cos%), (103) 

A"* A 

where Ap„ = Pn—Pn is the difference between the bubbling solution eigenvalue moments 
and the AdS^ x S^ moments (we have also encountered this difference of moments in the 
computation of the correlators in the gauge theory (17T1) ). In deriving these formulas, 
we have taken 

Po-pJi = 0, pi-p? = 0. (104) 

The first condition comes from fixing the radius of the two geometries, so that both 
are asymptotically AdS^ x S^ with the same radius of curvature. The second one is 
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the "center of mass" condition. One can also see from expanding Wigner's semi-circle 
law - which controls the eigenvalue distribution of the AdS^ x S^ vacuum solution - 
that p° = 1, p? = 0, p° = A/4, pO = 0, and p° = AV8. 

The prescription in [l5l HI] is to express the near boundary expansion of the metric 
and the other bulk fields in Fefferman- Graham form: 

dZ' dJ£dX^f (^..^2^ (^. 

+ Z^ (G'(4).,(X) + logZ2/i(4),^.(X)) + . . 

$2(X, Z) = Z' (logZ2<l>?o)(X) + $?o) W + •••), 

^\X, Z) = Z^-^<l>fo)(X) + . . . + Z'=<l>f2fc_4)(^) + • • • for fc > 2 , (105) 

where Z is the Fefferman-Graham radial coordinate and X* are coordinates on the 
boundary. The first terms in these equations, G(o)jj5 "^(o); and $/qn, are the sources for 
the stress tensor and the chiral primary operators of the field theory, while G(4)jj, $?qw 
and $^2/0-4) ^^® ^^^ normalizable modes of the fluctuations. For the bubbling solutions 
used in this paper, the non-normalizable modes that introduce sources vanish. 

Our first task is to perform a near boundary expansion of the bubbling supergravity 
solution. To do this, we introduce a radial coordinates -R, which is related to the 
coordinate x appearing in the bubbhng solution fl84|) by 



X = log (r + VWTl) . (106) 

In these new coordinates, the AdS^ x 5*^ metric (|96|) is given by 

ds'^ = 8c^ ( {R^ + l)dslas2 + -B^ — T + ^^^^2 + dy^ + cos^ y d^A . (107) 

In this coordinate system the conformal boundary is at i? ^ 00, where the metric on 
it is that of AdS2 x S'^. 

After some calculations, we get the asymptotic form of the deviations in fl99|) up to 
order O {R~^) terms: 

Ai = -— (4c4 + tanydyC4,) "H^ - 7T (Scs + tanydyC^) — 

~ 9fU8 t^^*^^ "^ ^^ ^^^ ^ ^y'^^ ~ ^^'^^ - 12 tan y dycl 

— Stanly {dyC4)'^ + 64c4(2cos2|/ — 7) + 32tan|/(cos2y + 2)dyC4\ — , 



i?4' 



1 11 1 

A2 = -— (4c4 + tan ydyCi) -^ - -^ (Scs + tan ydyC^) — 
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~T^7T7^ [96c6 + 16 tan y dyC^^ - AScl - 12 tan y dycj 

-3 tan^ y {dyCiY + 64:Ci{2 COS 2y + 1) + 32tan|/(cos2y - 2)(9j^C4] -^, 

^4 = — (4c4 + tanydyCi) -^ + -rr (Scg + tanydyC^ 



1 



+ 777-7^ [96C6 + 16tan?/9j/C6 - I6C4 -Atanydycl 
zU4o 



tan^ y{dyC4Y + 64c4 sec y cos 3y + 32 tan ?/ cos 2y dyCi\ 



1 11 1 

Ap = — (4c4 + tany 93/C4) ■^ + -^ (5c5 + tany SyCg) -^ 

~^9ru8 t^^*^*^ ~^ letanySj^ce - IQcj -Atanydycl 

— tan^ y{dyCi)'^ — 64c4(2 cos 2?/ + 5) — 32 tan y (cos 2?/ + 2) ^^04] — . 

(108) 

We recall that the functions Cm{y) are given in terms of the moments of the eigenvalue 
distribution by (11031) . Therefore, we have written the deviations in terms of the matrix 
model data, the various moments of the eigenvalue distribution. 

One finds similarly that the RR four-form deviation fl90|) is given up to order 
O (R-^) by 



= 1 + 7T (I6C4 + 8 tany 9j^C4 - dyCi) ;m + 7^ (^0^5 + 9 tanydyC^ - dyC^) — 

+ — — [24c6 + 10 tan y dyC^ - dlc^ - 96 sin^ y C4 
zoo 

+10sec?/(sin3|/ — siny) dyC^ + d'^c^il — secy cos 3y)] — . (109) 

Here we note that the coordinate R used in our expansions (I108|109l) is closely 
related to the Fefferman-Graham radial coordinate Z defined in (11051) . The relation is 
given by 

^=;^-^ + C?(l/^^). (110) 

This relation allows us to use compute the correlators for the various local operators 
up to dimension 4 by isolating the relevant term in the 1/ R asymptotic expansion 
of the corresponding bulk fluctuation, even though R is not the Fefferman-Graham 
coordinate. 
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Now that we have the exphcit form of the deviations that we need to calculate 
our correlators, we expand the deviations in a basis of spherical harmonics of S^. We 
decompose the metric and the RR 5-form into an AdS^ x 5*^ part and a perturbation 

gMN = 9mN + h-MN, 

F, = Ft^ + U (111) 

In this subsection, we are interested only in the fluctuations of the metric and of the 
5-form on S^, which only contain scalar harmonics |57)^^l 

f abode = 2_^ ^ ^ ^abcdeY , 

hab) = J2<pU^ia^b)Y' + ..., (112) 

where a, b, . . . are indices on S^, A^ = —k{k + A) is the mass of the k-th scalar spherical 
harmonic, and e is a volume form of S^ with unit radius. The brackets on the indices 
instruct us to take the symmetric and traceless combination. In the last line we omit 
terms involving vector and tensor spherical harmonics that are irrelevant for us. Later, 
for the computation of the stress tensor, we will also need 

K. = J2hlY', (113) 

where fi and z/ are indices on AdS^. 

Since the circular Wilson loop preserves an 5*0(5) subgroup of the R-symmetry 
group of A/" = 4 super Yang- Mills, so does the bubbling supergravity solution. There- 
fore, in performing the harmonic decomposition of the solution, only SO (5) invariant 
spherical harmonics contribute. These spherical harmonics depend only on the polar 
angle of the S^, which we identify with the coordinate y G [— 7r/2, 7r/2] of the base. 
The metric on the sphere is given by ds"^ = dy"^ + cos'^yds'^4, as can be seen from the 
AdS^ X 5*^ vacuum solution f[Ml) . The 5*0(5) invariant spherical harmonics are given 
in this coordinate system b}o 

Y'{y)=Afjcf{cosy), (114) 

where Cj (cos y) are Gegenbauer polynomials and the normalization factors are chosen 
as in (A. 2) of [41^ and in ([5]) 



^^For the correlators computed in this paper, the fluctuations in the other supergravity flelds do 
not contribute. 

""^^A brief review of 5*0(5) invariant harmonics and of Gegenbauer polynomials can be found, for 
example, in the Appendix of [37] . 
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This fixes the normahzation of the 5*0(5) invariant spherical harmonics to 



3J! 

In particular, the explicit normalization for the harmonics we will need is given by 

Afo = l, M = i ^^ = ^—, Ar3 = ^^, Af^ = ^—. (117) 

4' 2^30 8v^ 20^7 



For the reader's convenience, we list the explicit form of the first few harmonics: 

Y^ = MiA sin y , 

Y^ =7V2(-2 + 12sin2y), 

Y^ = A/'3(-12sin|/ + 32sin^l/), 

Y^ = ^fi{3-A8sm^y + 80sm'^y). (118) 

More details about the expansion in fluctuations and about general properties of 
spherical harmonics are given in [l5l HI] . 



4.3 Chiral primary operators 

To compute the one-point functions for chiral primary operators Oj we need appro- 
priate combinations of the trace of the metric and of the RR four-form fluctuations, 
which are mass eigenstate of the Laplacian on the sphere [58] 

The flrst expansion in (I112p can be inverted to flnd vr^ by using the orthogonality of 
the spherical harmonics 



^ = rn/2 — (120) 



/- V2 ^y ^» ^'^ ^°^^ y 

nildyiYf^ycos^y 

where, in our case, h'^ = 4A4 -|- Ap, which follows from (1991) . More precisely, we need 
to pick the appropriate terms in the expansion (11081) . The terms to compute depend 
on the dimension J of the dual chiral primary operator Oj. For O2 we have to select 
the coefficient of the R~^ term in (I108p . for O3 the coefficient of the term R~^ in (11081) . 
These are the same as the coefficients of Z^ and Z^, respectively, due to ( JllOp . The 
coefficient of the Z'^ term in the expansion of the quantity A is usually denoted with the 
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notation [A]k. In general [A]k differs from the coefficient of f? '^o From tlie explicit 
expressions in (11181) and from the expansion for A4 and Ap in (11081) we get 



M4 = ^ (-84AAP2 - 45(Ap2)' + IOOAP4) , (121) 

where hatted quantities in these formulas denote gauge invariant quantities at first 
order in the fluctuations. Note in particular that the scalars entering in the dimensions 
2 and 3 computations are automatically gauge invariant, whereas this is not the case 
for dimension 4 operators, as we will see presently. 

Similarly, we can invert the second equation in flll2p to get b'': 
^^ J2 - iH2 - J(y2)2eos4^ -^T^P' ' 



J[A pJ,{-j^)Y^cos^y _ 2VTO 
JiY^ycos^y 3A3/2' 



r,31 _ rTsi _ J i^Fy\2 {-YlJ ' ^Ub y _ ^ViU 

[0 J3 - [6 J3 - rfV3\2^.^^4.. - ~T^Jn^P3 ' 



Again, dimensions 2 and 3 quantities are already gauge invariant, unlike 6^. 

Gauge invariant combinations for the k = 4 case can be nonetheless easily formed 
using 0f w as explained in [l5l HI]. This coefficient can be obtained from the third 
expansion in flll2p . using some standard properties of spherical harmonics 

J D^'^D'^h^^.^Y' = 4)4 (1 + ^) A^ |(r^)^ (123) 

so that, integrating by parts and using fl99l 11081) we get 

^'^(«)J^ 4 (1 + f ) A4 /(F4)2 

1 / ([A,]4C^gy^ - 4[A4]4tanyg,y^ + f [A^ + 4A4]4y") 

4 (1 + f ) A4 jiY^r 

- -^. (124) 

One can now form the following gauge invariant combination of fluctuations [H] 

[^14 = M4 - ^' WJ4 = ^ (20^P4 - 9(Ap2)^ - 2OAAP2) , 



^^ For the quantities relevant to our particular computations, the differences cancel out. 

36 



Pl4 = [b\ - \ [ct>%^]^ = -^(AP4 - AAp2). (125) 

We have at this point all the ingredients to construct the gauge invariant mass eigen- 
functions to linear order in the fluctuations: 



M. = ^([*']^-«°i'^w = ir'^''^'' 



M4 = I^(P']4 - 80PI4) = ^, (20^ft - 3(AP2)' - 2OAAP2) . (126) 

Using holographic renormalization we are now able to extract the one-point func- 
tions of various local operators. Given the local five- dimensional supergravity action 
together with the associated counterterms, correlation functions can be computed by 
differentiating with respect to the non-normalizable modes(=sources) of the bulk fields. 
In the gauge theory, we have computed the correlation functions of unit normalized 
operators ([3]). The one-point functions of the unit normalized chiral primary operators 
O2 and O^ are given in terms of the supergravity fluctuations by [S] 

{02ix))^ = j^s% = y Y 3AP2, 

while the expectation value for the dimension 4 operator contains non-linear terms, as 
anticipated above, and reads 



(O^ix))^ — 



2 / 2\2 



2s^ + ^--a422 (s^) 
3z(4) ^''^ ^ 



27r2 5 
with the triple overlap function 0422 being 



-^(Ap4-AAp2), (128) 



a422 = ^ / Y\Yy = ^ (129) 



1 [Y\Y^f = ^ 
71^ J ^ ^ 800 

and z{4) = 1/240 as defined in 0117p . These are the final results from supergravity for 
the one-point functions of low dimension chiral primary operators. 

The non-trivial information about the correlator is in the function S^^^j defined in 

m 

{Ojix))^ = EnjYie). (130) 
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The dependence of the correlator on the choice of representative of the chiral primary 
multiplet factorizes, and it is captured by the spherical harmonic function Y{6), where 
6^ determines the coupling of the scalars to the loop ([1]). Since the bubbling super- 
gravity solution is S0{5) invariant, the supergravity correlator computes the one-point 
function of a chiral primary operator which is 5*0 (5) invariant and can be constructed 
from the S0{5) invariant spherical harmonics Y = Y^ in (11141) . For the choice of scalar 
coupling we have made where 6* = (1, 0, . . . , 0), we have that 



■.;,.,„^_.^(^ + 3)!_,/(•^ + 2)(•^ + 3) 



Therefore, we can compute H/j^j by dividing fll27|128l) by fll3ip and obtain 

Si?,2 = V2yAp2, S^,3 = 2^^^Ap3, S^,4 = 2^(Ap4 - AAp2). (132) 

The agreement between the supergravity results (I132p and the gauge theory computa- 
tions (1711) is then manifest! 



4.3.1 The small representation limit 

In comparing the gauge theory and supergravity results we have not used the explicit 
expressions for the moments p„. Here we evaluate them for the rank k symmetric 
and antisymmetric representation respectively and for the specific case of J = 2. For 
these small representations, the bulk computation of the correlator can be performed 
in terms of a probe D3 and D5-brane respectively [37]. Even though the curvature gets 
large in the interior of the bubbling geometry, the curvature is small near the boundary, 
which is sufficient to compute these correlators. 

From fll27p above and from fll72p and fll77p in Appendix O, we have (after trans- 
forming to M^) 



,2 



{WR{e,a)) L^WS^P' 



. TTA/l'^Vl + K? symmetric case , 
{Wnie,a)02iL)) Aa' N 5 ' ^ ^' 

- 1AP2 = ^ 

4a_ i^^sin^Ok antisymmetric case . 



L* Y 3 67r 

(133) 
On the other hand the probe D-brane computation in [37] gave 

{Wn{9,a)02{L )) a^ 
{WR{e,a)) 

with 



|lC.,.,y-(^). (134) 



9 79 \ 

cs,2 = 4V2WI + fi:2 , CA,2 = -^ sin^ Ok . (135) 

OTT 
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Putting everything together and recalhng that Y'^{tx /2) = a/5/6 one finds exact agree- 
ment with (I133p . 

4.4 Stress tensor 

We now move on to the computation of the one-point function of the stress tensor. 
As explained in [15], HI], we need to compute (l -|- |7r°) g'^^ + /i°^. Here /i°^ are the 
zero-modes of the metric perturbation on AdS^ (recall that Y^ = 1) 

8 r/^ 

hlu = IT ^^^M^ cos^y, (136) 

which explicitly read 



5^ \ A„2 

AdS2 



hld.^d.^ = [R^ + 1) ( '^"^g^gr""^ + 0{R~')) ds 



24A2i?4 ^ V i?2 + 1 

. /Ap2(45Ap2 + 14A) , ,^,^_5^\ . 2 



+R' { " ^ 2^^;^; — - + 0(^-') ) ^4^ , (137) 



while 7r° is the zero mode of the trace of the metric perturbation h^ 

"^2 9.c;i'An„^2 

t/2 

Then, the modified metric 



7r° 



— r^' t/y /i^ cos^ y = - ^^^ft^f + OiR-') . (138) 

3vry_w2 8A2i?4 ^ ^ ^ ' 



ds' = (^1 + ivr"^ <7j, + hlJx'-dx'' 
is expanded in 1/R as 



(139) 






12A2i?4 ^ V ;y 2 

. («^ . 1) (l . f) ,,. ,,,^ , (i , g) _|^ , «. (i , H) ,.. . (140) 

We now introduce the following near boundary coordinate 

1/^1 1 2 , P2 4 



«-^-4-+r*l- (1«' 
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so that the metric becomes 



•2 ■? I 16 4 



+ ^-^+h^ + ^+P3h=M<i4- (142) 



Note that z is the near boundary coordinate for the modified metric (1139p and thus 
differs form Fefferman- Graham Z introduced in flllUp . 

From [H], the stress tensor correlator is given bjo 

+ g[Tr^?2) - (Tr^(2))']^(o)^i " 2 ^2))^^ + 4^{2)yTr^(2) j , (143) 

where the g(k)ij^s are the analogues of the quantities in (11051) for the metric fll39p . 

Plugging fll42p into the expression above and using [s^]2 = ^^Ap2, which appears 
in the dimension 2 chiral primary calculation, we get 

{T,j{x))wdx'dx^ = ~^^Ap2ids\as2 - ds^) + ^^i.ds\dS2 + ^^lO- (144) 
This precisely agrees with the gauge theory computation ( |72l) 

N N^ 

including the conformal anomaly contribution ( TT2l) 

A^2 
{Tij{x))wdx'dx^ = hw{ds\^s2 - ^s^) + ^^(dsl^s^ + ds^) . (146) 
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""^^Note that </(„) denotes the coefficient of the z" term in the expansion of (|142p . i.e. 

ds^ ^ ^ + ^(5(o)-y + • • • + z"g^„^,j)dx'dx^. 
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A Weyl transforms between boundary metrics 

In this appendix we discuss the two Weyl transformations relating W^ and AdS2 x S"^, 
which we have used in section 12.11 The first transformation is relevant for the circular 
loop computation, while the second one for the straight line. 

Let us parametrize M^ using two sets of polar coordinates so that 

^4^ = dr^ + r^di/j^ + dL^ + L^d^l (147) 

These coordinates are relevant for a circular loop, which we take to be defined by r = a 
and L = 0. By making the following change of coordinates 

(r2 + L2-a2)2 + 4a2L2 a^ ^ . ^ r ~ ■ n 

r = -^5 = 7 7^ , r = r smh p , L = r sm fc' , 

Aa^ (cosh p — cos u)"^ 

(148) 
we find the metric 

dsl, = ^ [dp" + Sinh^ p d^2 + rf^2 ^ gi^2 Q ^2^ ^ ^^4g) 

which is conformal to ^4^5*2 x S*^ in global coordinates. Note that the conformal factor r 
is that in ^ and that the loop, which was located at r = a, L = in M"^, gets mapped 
to the conformal boundary of AdS2 x 5*^, namely the boundary of the Poincare disk. 

Now, under the conformal transformation (I149p a dimension J operator Oj trans- 
forms as follows: Oj — >■ r~^Oj. This proves the relation between the form of the 
correlator in M"^ ([8]) and the one in AdS2 x S"^ ([7]). 

The metric for M^ can also be written as 

^44 = dt^ + df" + l^ds\2. (150) 

We place the straight line at / = 0. In this case the Weyl transformation to 74^5*2 x 5*^ 
is simple: 

asjg4 = I "S^^^xs* ' (151) 

where 

dsAdS2xs^ = p + dss2 (152) 

involves the AdS2 metric in Poincare coordinates. The operators transform as Oj -^ 
l~^Oj when going from AdS2 x S"^ to M^, thus proving (TTOl) . 
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B Relating the stress tensor to a chiral primary via 
a GL twist 

Here we rederive the relation (13T]) between the correlator of the Wilson line with the 
stress tensor and the correlator of the Wilson line with the dimension two chiral primary 
O2 from a different point of view. 

The basic observation is that the supersymmetric Wilson line is closed with respect 
to the BRST charge of the Geometric Langlands (GL) twist |l9]o Since, by definition 
of a topological field theory, the twisted stress tensor T' is BRST exact, it follows that 
the expectation value of T' in the presence of the Wilson line is zero: 

{Wnie, a) T;,) = {WR{e, a) {Qgl, V,,}) = -{V,, {Qgl, Wr{9, a)}) = . (153) 

We can then compute the difference between the stress tensors of the twisted and 
untwisted theories and consider its correlator with Wr{6, a). This will turn out to give 
the wanted relation with the correlator of 02- 

We do not need to consider the kinetic term of the gauge field since this is an R- 
symmetry singlet not affected by the twist. We also ignore fermions at first. Before 
the twist the action in a generic curved background is 



9ym J V D / 



(154) 



where, as before, i = 1, . . . , 6. We now identify an 5*0(4) C 5*0(6) with the Lorentz 
group, so that i = fi,5, 6, and define a = 75(0^ + «0^). The GL twisted action is then 
given by (see equations (3.46 - 3.48) in [49j) 



9ym J ^ 



iJYM J ^ ^ 

+2D^(7D^a- 2[0^,a][0^a] + [a,af) . (155) 



The covariant derivatives D^ contain both the gauge and metric connections. In flat 
space S = S'. 

Let us now compute the stress tensor by taking the variation of the action with 
respect to the metric and setting in the end g^i^ = Sfj,^. One finds 



T,u ^M--^^(^ ^') 



Qfiiy — 0^1^ 



-"^^More precisely, the GL twists form a 1-parameter family of twists, where the parameter, t, is the 
projective coordinate on CP^. The Wilson line is closed with respect to the t = i twist. 
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2 „ /1„ „o, ,. ,.-. 1 



Tt(^-5,,D\<P'<J)') - -D^D,{<j)'<p') + 0(^DV.) + ^p(0(M^, 



9ym ^Q 6 



1 1 

-Dp{(f)pD(f,(f)^)) - -5^uDpD^{(f)p(f)^) - -D'^{(t)^(t)y) + DpD(^p{(l)^)(l)p) 



+ [<!>,,, <t>p][<t>uAp] + '^[<P{t,,^][<t>^),(^]) , (156) 

where we have used the following formulas 

K-^ = \9"'{Dp6g,, + DJgp, - DJgp,) , 
6R = Rpjg^'' + gp.D'Sg^'' - DpDJg^'' , 
SRpa = ^g,,D^pD^)6g^'' + ^gppg^,D^6g^'' - gp(pD,D^)6g^'' , (157) 

and integrated by parts. Imposing the equations of motion 

we can eliminate the quartic terms in (11561) and arrive at the final expression 

Tpu - T;, = ^ Tr (l:Sp,D'i<P^<P^) - \d pD ,{<P^ <P^) + Z},(0(^D,)0,) 
9ym V " " 

-Dp{(f)pD(^p(j)^)) - -Sp^DpD„{(f)p(t)„) - -D'^{(t)p(t)y) + D(^pDp{(j)^)(t)p) \ . 

(159) 

Note that on-shell d^{Tpy — T'^^) = 0, as it should be. 

Let us take now the Wilson loop to be a line along the /i = 1 direction and consider 
the correlation function with the stress tensor. By 5*0(5) symmetry we can say that 
(VrR(^,a)Tr(0>^)) = (Vrij(^,a)Tr(020^))(5"^ and (VrK(^, a)Tr(0V)) = where a = 
2, . . . , 6. Moreover, derivatives in the /^ = 1 direction vanish because of translational 
invariance and {Wii{d,a)TT{(j)°'dp(l)^)) is non zero only if a = /3. Putting everything 
together we have (here a = 2, 3, 4) 

{Wn{e, a)Tn) = -7^^9,9, {Wnid, a)Tr (0V' - V) > , {WR{e, a)T,,) = 

{Wnie, a)Tab) = -:^idadh - 6abdcd,) (Wnie, a)Tr (0V' - 0'0') > (160) 

'JOym 

and it is also immediate to realize that {Wr{9, a)T^) = 0. 

The operator Tr {(f)^(f)^ — (f^cf^) is a chiral primary of dimension two. The depen- 
dence (oc 1//^) of {W r{9 , a)Ti {(f)^ (f)^ — (f^cfP')) is consistent with the functional form 
(US]) of {WR{e,a)Tp^). One then finds that 

4 TV 

'Jfi'yM 3v2vr^ 
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Here we have used ¥{6) = j^{0^0^ - 0^6^) evaluated at ^ = (1, 0, . . . , 0). 

We see that the bosonic contributions have already reproduced the full result (ETl) . 
Thus it should be possible to show that the fermionic contributions sum up to zero, 
though we do not perform this computation here. 

C Moments in the small representation limit 

In this appendix we compute the explicit expressions for the second moments (^^) of the 
hermitian eigenvalue distributions in the cases of rank k symmetric and antisymmetric 
representations. Wilson loops transforming in these representations are described by 
D3 and D5 probe branes, respectively, having k units of string charge dissolved in their 
worldvolumes 



C.l Symmetric case 

Let us call the eigenvalues ^i and label them in increasing order, ,^i < . . . < .^at. Without 
the Wilson loop insertion, the eigenvalues are distributed on the interval [— vA, "\/A] 
and satisfy Wigner's semi-circle law derived from the saddle point equations 

- ^^i'^ + E ,(0) ! ,(0) = foralH = l,...,iV. (162) 

.(0) 

is given by c.o(C) = g'vM E. 1/(C - d°^) = 2C - 2 V'C^^A. 



As is well-known, the last eigenvalue is ^]y = vA and the resolvent of the matrix model 



Inserting in the path integral a Wilson loop in the rank k symmetric representation 
consists in moving the last eigenvalue ^at a distance k away from the interval. The 
saddle point equations become then 

47V 2 2 
-^^^+ V - + - - = for^ = l,...,iV-l, (163) 

^ ,y,,y^ ^^ - ^J ^^ - ^^ 

AM ^^-^ 9 

We make the ansatz that 

e. = d°^+56, SC, = 0{1/N) forz = l,...,iV-l, (165) 

while the shift for the ^n is large as remarked above. Then (11641) implies that 

^^ = y/x^lT^ + 0{1/N) , (166) 
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where n = yXk/AN. We want to compute 

7V-1 N 

= ^ E ^i°^^^^ + ^^^ - J^^n' + 0{l/N% (167) 

By taking the difference between the two saddle point equations fll62p and fll63p 
for z = 1, . . . , A^ — 1, we get 






We multiply this equation by Q and sum over i from 1 to A^ — 1. By playing with 
the dummy indices, one can show that 

N-l ^^ „^ Af-17V-l 



^^i - ^0 = _2 \^ \^ ^^' 

(0) _ A0)^2 ^ ^ t(0) _ t(0) 



9 \^ f (0) Y^ "^» ~ "^i _ _9 \^ \^ 

^Z^^i Z^ ,.(0) ^(0)no ~ ^2^2^ 



ItK^) _ ft.'JA2 -^^ -^^ ti>JJ _ t 

i=\ j^iJ^N \^i ^j J i=l j=l Si Sj 

N-l 



iiV^^(o) 



5^er5e. + 0(i). (169) 



i=l 

(0) 



In the second line we used the saddle point equation for Q ' . Thus 
8N v-^ ^rni^. v^ :.(n^ / 2 2 



- X E «.'°'*«. + E «.'"' -Twtjm + TWV I = 0( VA') . (170) 



i=l i=l \ Si SAT ?i <,N 

from which follows 

TV-l 



I 5^ ef^^e. = ^(-e^c.o(^^) + d°^^o(d°^)) + 0(l/iV^). (171) 

1=1 

By collecting everything, we get 

A 



(e) - {e)o = -^vTT^+o{i/N'). (172) 

C.2 Antisymmetric case 

The effect of inserting a Wilson loop in the antisymmetric representation is to create 
a hole in the [— vA, vA] interval so that the distribution splits into two groups with 
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k and N — k eigenvalueso The shift is 0{1/N) for all of them and the saddle point 
equations read 



-— e. + V- — - = o ioTt = i,...,N-k, 
^ .y. ^' ~ ^' 

^^-6 + 1 + 5Z ^4^ = ioTt = N-k + l,...,N. (173) 



A ^6-0 



Subtracting fll62p from these expressions one gets 

-X^e. + 2$:-j|f^ = 0(l/iV) for. = l,...,iV-fc, 



(174) 



—6^, + 1 + 2 ^ J^^_ J^l = O (l/N) ioTz = N-k + l,...,N. 



In this case the expression for the second moment reads 

9 ^ 

{e)-{e)o = j^j2^r^5^^+o{i/N^). (m) 

4 = 1 

Multiplying (11741) by ^^ , summing over the respective ranges of z, and finally summing 
the two equations one has 

N N c^ c^- N 



^d°'^6 + 2Ed°'E ,;g"X + ^ «f' = 0(l/iV). (176) 



Using (11691) and defining ^ = vAcos6', one has up to orders O {1/N'^ 



i=l i=Ar_fc+l ^>ik 

\3/2 /-efc X3/2 

= / desin^ecose = sln^^fe. (177) 



^°Note the difference between the symmetric and antisymmetric representation: the former can 
have arbitrary rank fc, consistently with the fact that we can move ^n arbitrarily far away from the 
interval, whereas the latter must have k < N , with the hole confined inside the interval. 
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